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CATALOG AND ATLAS FROM A SENSITIVE 1.5 GHz RADIO SURVEY 
AROUND THE NORTH ECLIPTIC POLE 


INTRODUCTION 


The correlation of deep multifrequency surveys is useful in studying the relative 
populations and cosmic evolution of various classes of extragalactic objects (e.g., active 
galactic nuclei, starburst galaxies) and other aspects of observational cosmology, such as 
galaxy clustering and the origin of the X-ray background. The North Ecliptic Pole (NEP), 
located at a = 18*00™ and 6 = 66°30’, is a region of special importance for such studies. All- 
sky surveys from satellites often scan great circles perpendicular to the Earth’s orbit in order 
to keep the solar panels facing the sun, resulting in sensitive coverage of the ecliptic poles. 
Our survey was motivated by the very deep coverage obtained by ROSAT (Trimper 1983; 
Voges 1992) but is also useful for comparison with other satellite based surveys. The NEP 
region has been previously studied with radio surveys using the Green Bank 91-m telescope 
at 1.4 GHz (Condon & Broderick 1985, 1986; White & Becker 1992) and at 5 GHz (Becker, 
White, & Edwards 1991; Gregory & Condon 1991) and the Effelsberg 100-m telescope at 
2.7 GHz (Loiseau et al. 1988). Very Large Array and Ryle Telescope observations have been 
made of sources selected from the 38 MHz 8C survey (Rees 1990; Lacy, Rawlings, & Warner 
1992). 


We present here radio observations obtained of the North Ecliptic Pole region with the 
NRAO! Very Large Array (VLA) radio telescope (hereafter the VLA-NEP survey) that are 
10 -— 100 times deeper than these previous surveys. This survey was undertaken with the 
specific goal of correlating faint radio sources with faint X-ray sources observed in the recent 
ROSAT all-sky survey (Voges 1993; Brinkmann 1993). We observed 114 fields covering 
29.3 square degrees within ~ 3 degrees of the NEP. We have cataloged 2435 sources, with 
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about 9% fainter than 1 mJy. The 114 fields do not have uniform sensitivities (fields closer 
to the NEP were observed longer to match the increased sensitivity of ROSAT), and the 
sensitivity varies considerably across each field, but no biases have been introduced. The 
resulting catalog should be representative of the faint radio source population about the 


NEP. 


OBSERVATIONS 


Observations of the NEP were made on 31 December 1990 and 2 January 1991 with the 
C-configuration of the VLA. Details of the observations and the data analysis are given in 
the Appendix, and we briefly summarize them here. 


Pointings were made of 113 fields around the NEP, each separated by 30’, plus one 
calibration field 3.6° from the NEP. Typically three short exposure “snapshots” of 8-minutes 
each were made for each of the 28 fields centered within 1.5° of the NEP, while two 4- 
minute observations were made of the 85 fields centered between 1.5° — 3.0° of the NEP. 
The nominal resolution of these observations is 20”. In total, 29.3 square degrees were 
surveyed. Figure 1 presents a schematic of the regions about the NEP covered by each field 
and shows the corresponding numerical designation adopted for each field. Table 1 presents 
the observational parameters for each of the 114 fields observed. The columns in Table 1 list: 
(1) the field designation; (2)-(3) the coordinates of the field center at epoch J2000; (4) Orms, 
the root-mean-square noise level at the field center in mJy; (5) the minimum signal-to-noise 
ratio required for source detection; and (6) the major axis, minor axis, and position angle, 
respectively, of the restoring beam. The elements in columns (4)-(6) are described below 
and in the Appendix. 


Calibration, mapping, and flux measurement operations were performed using tasks in 
NRAO’s Astronomical Image Processing System (AIPS). Fourier transforms of the data were 
made, and the resulting maps were averaged in the image plane. The root-mean-square noise 
of each field (rms) was measured over a central region, avoiding bright sources. Each field 
was searched for sources with peak fluxes exceeding 6 x Orms- In some fields the minimum 
flux density was increased to 7-8 x Opms due to systematic errors in the image. Peak and total 
flux densities were measured for every source which exceeded the minimum signal-to-noise 
ratio for that field. 


The flux densities measured from images produced by a synthesis array require 
corrections. We have corrected the measured peak flux densities for each source in the VLA- ; 
NEP catalog for (i) the primary beam response of the antenna elements and (ii) chromatic 
aberration due to the finite bandwidth, and we have corrected the measured total flux 
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Figure 1 — (Top) The location of fields imaged around the North Ecliptic Pole (NEP). The 
two darkened regions are areas that were not imaged because of mispointed fields. The field 


around the calibrator 1748+700 is included with the 113 fields centered within 3° of the 


NEP. (Bottom) The numerical designations for the 114 fields. 


densities just for the primary beam response. Uncertainties in the corrected peak and total 
flux densities were estimated from the quadratic sums of several terms. These terms included 
(i) statistical error from the root-mean-square noise level, (ii) uncertainties in the calibration 
and flux measurement process, (iii) pointing errors in the antenna elements, (iv) uncertainty 
in the primary beam correction, (v) uncertanty in the bandwidth smearing correction, 
and (vi) uncertainties due to map artifacts. All of these corrections and uncertainties are 
described in detail in the Appendix. 


The positions measured from synthesis array images likewise require corrections for 
instrumental effects. The raw map must be expanded by a small factor to recover true 
positions. Positional uncertainties were taken to be the quadratic sum of two terms: 
(1) statistical uncertainty from the finite signal-to-noise ratio of the source, and (ii) systematic 
residual uncertainty from the application of the radial expansion. These corrections and 
uncertainties are also discussed in the Appendix. 


THE VLA-NEP CATALOG 


Table 2 gives the resulting VLA-NEP source catalog. The columns list: (1) the source 
name at epoch J2000; (2) the field designation where the source is best detected; (3)-(8) the 
corrected source position at epoch J2000; (9) A(position), the positional uncertaintity in 
arcseconds; (10) p, the distance of the source from the field center in arcminutes; (11) Stota, 
the corrected integrated flux density in mJy; (12) AStoa, the error in the integrated flux 
density; (13) Speak, the corrected peak flux density in mJy; (14) ASpeax, the error in the 


peak flux density; and (15) sh/ Orms, the signal-to-noise of the detection. Fully resolved 
double sources with bright lobes have been included as two catalog entries, with the total 
flux being that for the entire object. 


The last column of Table 2 contains additional notes where (ID) denotes possible 
identification with an object from another survey (see Table 3 and below), (S) indicates 
a source which is suspect because it lies on an image processing artifact, (D) is a double 
source, either a blend of two comparable sources or two resolved peaks with a common 
envelope, and (E) is an extended source, either an elongated source or an object with 
an asymmetric or complex shape. In addition (D?) and (E?) indicate possible doubles 
and extended sources whose features appear at less than 8 X Orms. Also included is the 
angular separation (in arcminutes) between the peaks of double sources, the angular size (in 
arcminutes) of extended sources measured from the peak to the most distant feature minus a 
slight resolution correction, and the position angle (in degrees) measured from the brightest 
peak. We identified 124 double sources and 174 extended objects; they range in size from 
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0.1’ - 2.2’ with 70% of the doubles and 60% of the extended sources larger than 0.5’. There 
are 45 double sources where we have have both peaks listed individually in the catalog; 
in these cases the second peak is indicated in the last coulumn of Table 2 by the symbol 
“>>”. Contour images of 54 sources with interesting structure are included in the Appendix 
(Figure A4), and these sources are indicated in Table 2 with the (F) identifier. With the 
exception of 11 sources, the position of each source is indicated in the atlas (Figure 2); the 
sources missing from the atlas are indicated in Table 2 with (X). 


The VLA-NEP source catalog was compared with four previous radio surveys which 
cover all or part of the VLA-NEP survey area. These include the 1.4 GHz Green Bank 
survey (White & Becker 1992), the 5 GHz Green Bark survey (Becker et al. 1991), the 1.5 
and 5 GHz VLA and Ryle Telescope observations of the 38 GHz 8C survey (Lacy et al. 
1992), and the 2.7 GHz Effelsberg survey (Loiseau et al. 1988). The details of determining 
positional coincidence are given in the Appendix. VLA-NEP counterparts were found for all 
of the 1.4 GHz Green Bank sources, and approximately 90% of the sources from the other 
three surveys. The comparison with the 5 GHz Green Bank and 8C catalogs found relatively 
few with more than one VLA-NEP counterpart within the search radius, however 36% of 
the 2.7 GHz Effelsberg and 78% of the 1.4 GHz Greenbank sources have multiple VLA-NEP 
counterparts. Where the identification of a VLA-NEP source is not likely to be subject to 


confusion (> 80% of the flux arising from a single VLA-NEP object), we have determined 
spectral indices between the VLA-NEP flux densities and the 2.7 GHz and 5 GHz catalogs. 


In addition to radio surveys, the VLA-NEP catalog was compared to catalogs 
of extragalactic objects listed in the NASA Extragalactic Database (NED). We found 
42 galaxies (most brighter than 15 magnitudes) and 27 [RAS sources with VLA-NEP 
counterparts. VLA-NEP counterparts were also found for the galaxy cluster Abell 2280, 
the planetary nebula NGC 6543, and the supernova SN 1989P, although the latter is most 
likely a chance coincidence from the large search radius adopted. 


The results of these tentative source identifications are listed in Table 3. For each source 
listed in Table 2 with an “ID” note the columns show: (1) the VLA-NEP source name; (2) the 
VLA-NEP corrected total flux density in mJy; (3) the name of any 5 GHz counterpart from 
Becker et al. (1991) along with the cataloged flux density in mJy; (4) the name of any 2.7 GHz 
Effelsberg counterpart from Loiseau et al. (1992) along with the cataloged flux density in 
mJy; (5) the name of any 1.4 GHz counterpart from White & Becker (1992) along with the 
cataloged flux density in mJy; (6) the name of any 8C counterpart from Lacy et al. (1992); 
(7) the name of any NED counterpart that is not from the four previous radio surveys with a 
visual magnitude, if available from NED; (8) the spectral index between 1.5 GHz and 5 GHz; 
(9) the spectral index between 1.5 GHz and 2.7 GHz; and (10) a reference for sources with 


information from NED. 


THE VLA-NEP ATLAS 


An atlas of contour maps of the 114 fields included in the VLA-NEP survey, and listed 
in Table 1, are presented in Figure 2. Each contour map shows the central 42.7! x 42.7' 
region of the field. The image is a 512 x 512 pixel array, with each pixel being 5” square. 
Contour levels are scaled to the root-mean-square noise O;m, of the field. Contours are set at 
+£3.5, 6, 18, 54, ... times opms, where the progression is logarithmic with succeeding contour 
levels at 3 times the flux of the previous level. 


All sources from the VLA-NEP catalog (Table 2) which fall within a field are indicated 
on the contour map by a cross. A large cross indicates a source whose position was taken from 
this field. Small crosses are sources whose positions are taken from neighboring, overlapping 
fields. These are the raw, uncorrected positions. The offset of small crosses from contour 
centers is an indication of the difference in uncorrected positions from one field to another. 
Accurate positions should be taken from the VLA-NEP catalog (Table 2). 


DATA PRODUCTS 


Tables 1, 2, and 3 (total size 0.5 MByte) are available electronically from R. I. Kollgaard 
(rik@astro.psu.edu) or P. Hertz (hertz@xip.nri.navy.mil). Other means of transport can 
also be arranged. Also available is a FORTRAN program giving the flux sensitivity at any 
location in the survey region. This gives upper limits to any radio source, taking into account 
the variations in sensitivity due to different exposure times and primary beam attenuation. 
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TABLE 1 
VLA-NEP Survey Field Parameters 


Orms SNR Beam 

(J2000) (mJy) ("7 " *) 

0 663000 0.141 6.0 4925 -53 
55.0 670000 0.128 6.0 4926 -51 
58.0 660000 0.068 8.0 2518 -60 
55.0 673000 0.134 6.0 4928 -50 
55.0 663000 0.083 6.0 4925 -52 
55.0 653000 0.112 8.0 2217 -69 
10.0 680000 0.105 6.0 5813-0 
02.0 670000 0.092 7.0 4924 -53 
53.0 660000 0.054 6.0 2318 -58 
08.0 683000 0.072 6.0 4216 -4 
42.0 650000 0.062 6.0 2217 -69 
08.0 673000 0.199 6.0 §237 -47 
56.0 663000 0.103 8.0 2518 -59 
31.0 690000 0.072 6.0 3417 3 
44.0 653000 0.063 8.0 2217 -68 
24.0 643000 0.051 6.0 1917 -77 
30.0 680000 0.092 6.0 6514 4 
09.0 670000 0.099 9.0 2319 -56 
48.0 660000 0.044 7.0 2518 -61 
26.0 650000 0.054 6.0 2217 -70 
36.0 683000 0.067 6.0 5114 7 
01.0 640000 0.059 7.0 1916 -78 
21.0 673000 0.187 7.0 3431 -45 
$7.0 663000 0.086 75 2818 -21 
33.0 653000 0.183 7.0 2217 -69 
03.0 643000 0.087 6.0 2215 -25 
05.0 690000 0.101 7.7 3018 2 
50.0 680000 0.052 7.0 1414 -43 
16.0 670000 0.226 6.0 3218 -13 
43.0 660000 0.080 6.0 2417 -27 
10.0 650000 0.066 6.0 2415 -22 
32.8 700551 0.070 6.0 2119 -36 
35.0 640000 0.050 6.0 1916 -86 
04.0 683000 0.112 6.0 4315 -6 
34.0 673000 0.058 6.0 3018 -16 
58.0 663000 0.073 7.5 2518 -24 
0 653000 0.076 6.0 2517 -23 

0 643000 0.081 6.0 2516 -21 

0 690000 0.070 7.0 2719 -13 
10.0 680000 0.079 6.0 1514 -41 
23.0 670000 0.067 6.0 1716 -36 
38.0 660000 0.078 6.0 2220 -52 
54.0 650000 0.068 6.0 2115 -25 
09.0 640000 0.047 6.0 1916 -2 
32.0 683000 0.132 6.0 5015 -1 
47.0 673000 0.068 7.0 1616 -42 
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30 19 


15 14 
20 17 
22 18 
30 17 


28 17 
15 14 


16 14 
17 15 
32 18 
37 18 


50.0 680000 0.094 6.0 3117 -15 
F100 18 1916.0 653000 0.097 6.0 1514 O 
F101 18 1929.0 690000 0.161 8.0 2917 -33 
F102 18 2004.0 663000 0.069 6.0 2414 -22 
F103 18 2052.0 673000 0.119 6.0 3517 -10 
F104 182118.0 650000 0.093 Lo 1614 -1 
F105 18 2152.0 683000 0.106 6.0 2622 -18 
F106 18 2208.0 660000 0.101 6.0 1514 -40 
F107 18 2305.0 670000 0.091 6.0 4017 -5 
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F109 18 2410.0 680000 0.129 7.0 2420 -30 
F110 18 2505.0 663000 0.063 6.0 2416 -22 
F111 18 2605.0 673000 0.103 6.0 2719 -9 
F112 18 2703.0 660000 0.153 7.0 2315 -22 
F113 18 2812.0 670000 0.097 6.0 2919-2 
F114 18 3006.0 663000 0.052 6.0 2415 -21 


NOTES. -- Amplitude self calibration was used in the following fields: F006, 
F029, F032, F037, F049, F0S7, F0S8, F065, F070, F086, F089, F097, F101, F105, 
F107, F109. F018: image shifted to include bright source, 2.3” on western edge 
not searched; F041: image shifted, 4.6’ on eastern side not searched; FOS: 
image shifted, 1.8° on northern edge not searched, Orms = 0.135 mJy within 15° 
of the bright source VLA-NEP 1758.6+6637; F058: Opms = 0.121 mJy within 10° 
of VLA-NEP 1758.6+6637. 
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Figure 2 — Atlas of VLA-NEP Fields (continued). 
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A SENSITIVE 1.5 GHz RADIO SURVEY 
AROUND THE NORTH ECLIPTIC POLE 


Ronald 1. Kollgaard, Wolfgang Brinkmann, Margaret McMath Chester, Eric D. Feigelson, 
Paul Hertz, Patricia Reich, & Richard Wielebinski 


To appear in The Astrophysical Journal Supplement Series 


Abstract 


A 29.3 square degree region surrounding the North Ecliptic Pole (NEP; a = 18'00™, 
6 = +66°30') was mapped with the Very Large Array at 1.5 GHz to support the deepest 
portion of the ROSAT all-sky soft X-ray survey. The resulting VLA-NEP survey catalog 
contains 2435 radio sources with flux densities in the range of 0.3 - 1000 mJy, including over 
200 fainter than 1 mJy. The 28 fields of the inner 1.5° have noise levels o ~ 60uJy, and the 
85 fields centered between 1.5° - 3.0° from the NEP have o ~ 120uJy. The typical spatial 
resolution is 20” HPBW, and most positions are accurate to < 2”. Approximately 6% of 
the sources are found to be extended with size > 30”. We have compared the VLA-NEP 
catalog with four other radio catalogs made at lower resolution, as well as with the NASA 
Extragalactic Database, and find counterparts for ~ 18% of the VLA-NEP objects. The 
normalized, differential radio source count is in agreement with previous studies. Between 1 
-— 150 mJy the slope of the logN-logS relation is 0.68 + 0.03. 


Al Introduction 


The correlation of deep multifrequency surveys is useful in studying the relative 
populations and cosmic evolution of various classes of extragalactic objects (e.g., active 
galactic nuclei, starburst galaxies) and other aspects of observational cosmology, such as 
galaxy clustering and the origin of the X-ray background. The North Ecliptic Pole (NEP), 
located at a = 18"00™ and 6 = 66°30’, is a region of special importance for such studies. 
All-sky surveys from satellites often scan great circles perpendicular to the Earth’s orbit in 
order to keep the solar panels facing the sun, resulting in sensitive coverage of the ecliptic 
poles. Our survey was motivated by the very deep coverage obtained by ROSAT (Trimper 
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1983; Voges 1992) but is also useful for comparison with other satellite based surveys (e.g., 
HEAO-1, Wood et al. 1984; JRAS, Hacking & Houck 1987). The NEP region has been 
previously studied with radio surveys using the Green Bank 91-m telescope at 1.4 GHz 
(Condon & Broderick 1985, 1986; White & Becker 1992) and at 5 GHz (Becker, White, & 
Edwards 1991; Gregory & Condon 1991), the Effelsberg 100-m telescope at 2.7 GHz (Loiseau 
et al. 1988); and Very Large Array and Ryle Telescope observations of sources selected from 
the 38 MHz 8C survey (Rees 1990; Lacy, Rawlings, & Warner 1992). 


We present here new radio observations (hereafter the VLA-NEP survey) of the region 
that are 10 - 100 times deeper than these previous surveys, with the specific goal of correlating 
faint radio sources with faint X-ray sources observed in the recent ROSAT all-sky survey 
(Voges 1993; Brinkmann 1993). Our observations were made with the NRAO Very Large 
Array. We have cataloged 2435 sources, with about 9% fainter than 1 mJy. The 114 fields 
do not have uniform sensitivities (fields closer to the NEP were observed longer to match the 
increased sensitivity of ROS AT), and the sensitivity varies considerably across each field. 
But no biases have been introduced and the resulting catalog should be representative of the 
faint radio source population about the NEP. 


Details of the observations are presented in §A2, with a discussion of the image 
processing, correction for instrumental effects, and error determination. In §A3 we present 
the VLA-NEP catalog, which is compared with other catalogs in §A4, and in §A5 we 
determine the log N-log S radio source count distribution. 


A2 Observations 


Observations of the NEP were made on 31 December 1990 and 2 January 1991 with 
the C-configuration of the VLA. The data were recorded with dual 46 MHz bandwidth IFs 
(i-e., using the standard “50” MHz bandwidth) with an effective center frequency of 1489.9 
MHz. Pointings were made of 113 fields around the NEP, each separated by 30’. Typically 
three short exposure “snapshots” of 8-minutes each were made for each of the 28 fields 
centered within 1.5° of the NEP, while two 4-minute observations were made of the 85 fields 
centered between 1.5° — 3.0° of the NEP. The nominal resolution of these observations is 
20”. However, during the first day’s observations some data was lost due to a failure of 
the entire VLA. This eliminated one scan from several of the fields, with the result that 
ten outer fields have relatively high noise levels and elongated beams. Also one inner and 
one outer field were mispointed 20’ to the north, with the result that there are two small 
regions where no reliable source detections could be made. In total, 29.3 square degrees were 
surveyed. Figure 1 presents a schematic of the regions about the NEP covered by each field, 
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and shows the corresponding numerical designation adopted for each field. 


The visibility data were calibrated using the tasks supplied in NRAO’s Astronomical 
Image Processing System (AIPS), treating each day separately. The total flux density was 
tied to observations of 3C 286 and the scale of Baars et al. (1977). The source 1748+700 was 
used as the local phase calibrator. Individual fields were observed during only one of the 
two days, although the flux density of 1748+700 on each day differs by only 0.6%, indicating 
that the absolute flux scale is consistent between the two sessions. As 1748+700 is only 3.5° 
from the NEP, the field surrounding it has also been included in the catalog. 


A2.1 Imaging and Source Detection 


Fourier transforms of the data were made using the AIPS task UVMAP for the two 
IFs separately, and the resulting maps were averaged in the image plane. The maps are 
1024 x 1024 with 5” pixels. “Natural” weighting was used, giving the highest sensitivity 
with some loss of resolution. Deconvolution of the inner 512 x 512 pixels (42.7' x 42.7’) was 
accomplished with APCLN. The data were self-calibrated twice using ASCAL correcting only 
for the phases. Following self-calibration, the resulting images from each IF were averaged. In 
a few instances, the presence of a bright source (S > 100 mJy) required a further application 
of ASCAL with a correction for phases and amplitudes. However, since this causes a slight 
redistribution of flux from the fainter to the brighter sources, amplitude corrections were 
limited to a few fields (see notes in Table 1). 


The root-mean-square noise of each field (7) was measured over a central region chosen 
to avoid bright sources. Typically, the deeper fields within 1.5’ of the NEP have noise levels 
o around 60pJy, while the fields in the outer annulus, 1.5’ - 3.0’ from the NEP, typically 
have o + 120uJy. Noise levels along the edges and corners of each field are 3 — 4 times higher 
than at the field centers due to primary beam attenutuation and bandwidth smearing. 


The 42.7’ x 42.7’ region CLEANed was searched for sources at a level 6 x o in most 
fields. In some fields, a visual inspection showed systematic errors on the image and the 
minimum flux density for a detection was increased to 7 - 8 x o. Determination of the peak 
flux density and position of each source was made with the AIPS task MAXFIT, which fits a 
quadratic function to the 3 x 3 pixel area around each peak. Total flux densities were found 
with the AIPS tasks TVSTAT and IMSTAT which determine the sum of pixels in a small 
region around the peak. Regions were selected which, by visual inspection, best enclosed 
the source. In a few cases image artifacts result in large errors in the total fluxes for some 


of the faintest sources (< 1 mJy). 


Many regions > 15’ from field centers were included on more than one field. Where 
sources were detected on these overlaps, only those from the field with the lowest o were kept. 
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This eliminates the possibility of a single object being counted twice due to its appearance 
on more than one field. Figure Al is a schematic representation of all sources found over the 
entire survey region. 
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Figure Al — The distribution of VLA-NEP scurces found about the North Ecliptic Pole. 


A2.2 Flux Density Corrections 


The flux densities measured from images produced by a synthesis array must be 
corrected for: (1) primary beam response of the antenna elements; (2) time-average smearing 
due to the finite integration time; and (3) chromatic aberration due to the finite bandwidth 
(Bridle & Schwab 1989; Cotton 1989). Very roughly, these effects increase quadratically with 
distance from the image center. 


The primary beam expansion presented by Bridle (1989) has been used, where the power 
pattern of each antenna element, A(p), is given by 


1 2 2n 
—_— = Gn (9% 1 
A(p) »y ( 0) ’ ( ) 
with 
@ = 9.920378 x 107! 
a, = 9.956885 x 1074 
@2 = 3.814573 x 107° 
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a3 —5.311695 x 107!° 
a, = 3.980963 x 107??, 


where p is the angular distance from the field center (in arcminutes), and the effective center 
frequency vp = 1.4899 GHz. 


We have not corrected for distortions due to time-average smearing because the effect 
is always very small. The data were collected with an integration time of 30 seconds which, 
for our freqency and array combination, leads to reductions in the peak flux of only 1% 
for a point source 30’ from the field center (Cotton 1989). We have also not corrected for 
“CLEAN bias”, a small systematic reduction in the flux density of sources that is thought 
to occur when large fields of view are CLEANed from short observations (Condon, private 
communication; Cornwell, private communication). Because most of our fields were observed 
more than once, there should be little effect from CLEAN bias. In the worst case, for those 
few outer fields where we obtained only a single snapshot, CLEAN bias may have reduced 
the flux densities by ~ 2¢, although for our maps other systematic errors on these fields 
likely dominate this effect. 


For this experiment, where natural weighting was used, the effects of bandwidth 
smearing correspond to the “square bandpass, circular Gaussian tapering” case discussed 
by Bridle & Schwab (1989). The reduction in the peak flux of a point source is R(p) = 
ue / Speak, Where Speak and Soeak are the true and measured peak flux densities, respectively, 
and 


vr (7) 
R(p) = — erf | —], 2 
(0) =F ent (4 (2) 
with 
p= ef , 
Yo CUHPBW 
y = 2vin2, 


where the bandwidth Av = 46 MHz, @uppw is the resolution (in arcminutes), and erf() is 
the usual error function. For this experiment the effective bandwidth and frequency were 
constant, and the effective resolution varies, sometimes considerably, from field to field. We 
set Ouppw to be the average of the major and minor axes of the elliptical restoring beam 


used by the AIPS task APCLN. 


Therefore, the flux densities she and shia measured on the CLEANed maps have been 
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corrected as 


S ae es 3 
pee = Rip) Alp) (8) 
t 
Scot 
Stotal = Alp) (4) 


where Speak and Sota are the corrected source values. 


The sensitivity of the survey is illustrated in Figure A2, which shows the corrected peak 
flux densities as a function of distance from the field centers. 
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Figure A2 — The peak flux density (Speak) of cataloged radio sources, corrected for the 
effects of the primary beam and bandwidth smearing, as a function of p distance from the 
field center. 


A2.3 Flux Density Uncertainties 


Uncertainties in the flux densities (in mJy) were estimated from the quadratic sum of 
several terms (e.g., Condon, Condon, & Hazard 1982; André, Montmerle, & Feigelson 1987). 
The first term is the statistical error from the rms noise level, 


bo 
A(p)’ (5) 


O71 


154 


where Orms is the noise level of the map, p is the distance of the source from the field center (in 
arcminutes), A(p) is the primary beam response given by equation (1), and bis an empirical 
term. Following Condon et al. (1982) we have taken 6 = 1.52. 


The second term arises from uncertainties in the calibration and flux measurment 
process, which we take to be 3%, so that 


o2 = 3x 107? Speak, (6) 
where Speak is the corrected peak flux density given by equation (3). 


For the third term, which incorporates pointing errors in the antenna elements, we have 
adopted the functional form of Condon et al. (1982) but have taken the pointing accuracy 
of the antennas to be «€ = 0.1” (Perley 1992), so that 


o3 = 5.5.x 107* pSpeak- (7) 

The fourth term represents uncertainty in the primary beam correction and we have 
adopted the expression used by Condon et al. (1982) 

Oo, = 11074 p? Speak. (8) 

Following André et al. (1987) we include an error term for uncertanties in the the 

bandwidth smearing correction. The peak and total flux densities of 60 unresolved sources 

selected from the 28 inner fields were inspected and compared to the measured ratio of the 


peak-to-total flux densities with that expected from equation (2). We found no dependence 
on p and from the observed scatter adopt 


Os = 5x 10-7 Spe. (9) 


To account for additional uncertanties in the determination of the total fluxes for faint 
sources due to map artifacts and similiar effects, we have included an additional error term 


0 Stota > ImnJy 
* 10 
- 0.15tote1 Stotaa < ImJy. (10) 


The total error in the corrected flux densities is therefore, 


1/2 
ASpeak = (o3 +03 + 03 +073 +02) é ; (11) 


1/2 
AStota = (0? topt+oft+oe+ ey) ; ; (12) 
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where the primed values of 72, 73, and a4 denote the use of the total flux density instead of 
the peak flux density in the evaluation of these error terms. 


A2.4 Positional Corrections 


The positions measured from images produced by a synthesis array also need to be 
corrected for instrumental effects. The raw map must be expanded by a small factor to 
recover true source positions. As discussed by Mitchell & Condon (1985) these corrections are 
related to the finite bandwidth and arise from: inaccuracies in the true effective frequency, 
the frequency dependence of the primary beam, and the variation in flux density across 
the bandpass for sources with spectral indices a # 0. All of these effects increase with the 
angular distance of the source from the image center. We have adopted the overall correction 
used by Mitchell & Condon (1985) who observed with the same frequency and bandwidth 
combination as we did. The measured position of each source has therefore been expanded 
radially outward from the center of the image by Ap (arcminutes), where 


Ap = (3.8 x 1073)p + (3.3 x 107®)p°, (13) 
and p is the angular distance of the source from the field center (arcminutes). 


In order to confirm the validity of equation (13) compensating for radial positional 
distortions, we made two tests: (1) an external test with a comparison of VLA-NEP survey 
positions with those of other samples with high positional accuracy; and (2) an internal test 
of VLA-NEP survey positions for sources which appear on more than one field. For our 
external test we looked at two samples. Patnaik et al. (1992) have imaged compact radio 
sources with a positional accuracy of 0.012”. Four sources from their list were detected in 
the VLA-NEP survey. In addition, Lacy et al. (1992) observed 57 sources around the NEP 
with arcsecond resolution. However, as the Lacy et al. (1992) sources were chosen from the 
8C catalog selected at 38 MHz, they are primarily lobe dominated objects. We used only 
eight 8C sources with compact structures (< 5”) where the positions are likely independent 
of wavelength and beam size. 


The position correction was further tested by comparing repeated measurements in the 
“unmerged” list of VLA-NEP sources of ~ 1000 source positions in overlapping regions of 
adjacent fields. Sources with (uncorrected) positions within 15” of each other in adjacent 
fields were considered to be multiple observations of the same object. The average separation 
between these multiple entries is minimized when the survey positions are corrected by 
equation (13). 


Figure A3 shows the offsets between VLA-NEP positions and positions given by Patnaik 
et al. (1992) and Lacy et al. (1992) for twelve sources in their surveys as a function of the 
VLA-NEP source’s angular distance p from the field center. The curve shows the correction 
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Figure A3 — A comparison of NEP source positions, as taken directly from the images, with 
sources with very accurate positions from Patnaik et al. (1992, crosses) and the 8C survey 
Lacy et al. (1992, triangles). The large error bars represent offsets of NEP survey sources 
found on more than one field. The solid line is the positional offset predicted by the radial 
expansion of Mitchell & Condon (1985). 


given by equation (13). It is clear that the radial correction follows the observed offsets, 
though the scatter becomes large for p > 15’. Also plotted are the mean and standard 
deviation of offsets between multiple observations within the VLA-NEP survey in three 
ranges of p. These are also consistent with the expectations of equation (13), although the 
scatter is larger since positional errors are introduced twice from each of the two overlapping 


NEP fields. 


Additional non-radial distortions, such as those found by Taylor & Seaquist (1985), can 
arise from the use of a finite integration time and the neglect of the w-term (Cotton 1989). 
As mentioned above, the effect of time smearing is negligible for this array and frequency 
combination. The neglect of the w-term is more problematic. For the case of a single 
snapshot the use of a 2-dimensional Fourier transform (i.e., UVMAP or MX in AIPS) which 
cannot account for the w-term will lead to an apparent position shift, Ap, (arcminutes), in 
the direction of the zenith, and with a magnitude 


Ape = 1.5.x 10‘ p? tanZ (14) 


where p is the angular distance of the source from the phase center of the map (arcminutes), 
and Z the zenith angle of the observation. 
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We made three tests of the magnitude of distortions arising from neglect of the w- 
term, using fields which contained bright sources with accurate positions from Patnaik et 
al. (1992) and Lacy et al. (1992). First, the data were reimaged treating each short scan 
of a field independently. The positional offset of the test sources from each snapshot was 
compared with that expected by equation (14). Second, the data were reimaged using the 
AIPS task MX with a phase shift to move each test source to the center of the image. As the 
w-term distortion depends on distance from the phase center, this procedure will eliminate 
it for the test source. Third, the individual snapshots from each field were reimaged with the 
new AIPS task UVADC. UVADC takes clean components which have been subtracted from 
a uv data set, corrects their position for the effects of the w-term, and then adds the Fourier 
transform of the shifted clean components tack into the data. After each test, the radial 
correction in equation (13) was applied. Although application of the these methods did move 
the sources by 1-2”, the resulting positions were not closer to those of Patnaik et al. (1992) 
and Lacy et al. (1992). Thus, our positional accuracy does not appear to be limited by the 
w-term and we have not applied any w-term correction to the VLA-NEP survey positions. 


A2.5 Positional Uncertainties 


The net positional uncertainty Apos = /(Aa)? + (Ad)? was taken to be the quadratic 
suin of two terms. The first term arises from the use of the parabolic fitting routine MAXFIT 
and depends on the signal-to-noise ratio of the source and the size of the beam. Using the 
AIPS tasks IMMOD and CONVL we created simulated point sources with noise which were 
convolved to various beam sizes. MAXFIT was then used and the difference between the 
measured and true position found. In the range of beam parameters applicable here we 
emprically estimate the error due to MAXFIT to be 


_ 3] (25+ Oupsw o 
logo; = 3 ( 56 — log c : (15) 
pe 


where o is the measured noise of the map center, Oyppw is the effective resolution in 


arcminutes, and S. y is the uncorrected peak flux density. 


peak 


The second error term arises residual errors remaining after the application of the radial 
expansion correction discussed above, arising from instrumental or deconvolution errors that 
are not well understood. For sources with p < 15’ appearing in Patnaik et al. (1992) and 
Lacy et al. (1992) shown in Figure A3, the standard deviation of the net positional offset is 
1.2". For sources with p > 15’, we used the internal VLA survey comparison of objects in 
the “unmerged” list shown in Figure A3, giving a lo deviation of 2.7”. Since errors in two 
fields contribute to this deviation, the uncertainty of a single position is 2.7”/./2 or 1.9". 
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Thus, the second positional error term is 


1.2” p< 15’ 
— i ; 16 
ue 1.9" p> 15! (16) 
and the total positional uncertainty is 
1/2 
Dow = (02 bas). (17) 


AS The VLA-NEP Catalog 


Table 1 presents the observational parameters for each of the 114 fields observed. The 
columns list: (1) the field designation; (2) - (3) the coordinates of the field center (J2000); 
(4) o, the rms noise level at the field center in mJy; (5) the minimum signal-to-noise ratio 
required for source detection; (6) the major axis, minor axis and position angle, respectively, 
of the CLEAN restoring beam; and (7) a flag to denote fields where amplitude self-calibration 
was used. Brief additional notes have also been included at the end of the table. 


Table 2 gives the resulting VLA-NEP source catalog. (See §6 for details on how to 
obtain the entire catalog, we include only a sample page here.) The columns list: (1) the 
source name (J2000); (2) the field designation where the source is best detected; (3) - (8) the 
corrected source position (J2000); (9) Apo, the net positional uncertaintity in arcseconds; 
(10) p, the distance of the source from the field center in arcminutes; (11) Stora, the corrected 
integrated flux density in mJy; (12) AS,ota1, the rms error in the integrated flux density; (13) 
Speak, the corrected peak flux density in mJy; (14) ASpesx, the rms error in the peak flux 


density; and (15) she /o, the signal-to-noise of the detection. Fully resolved double sources 
with bright lobes have been included as two catalog entries, with the total flux being that 
for the entire object. (Slight differences exist between the total flux values for most doubles 
because of variations in the flux density corrections due to the different distances of the lobes 
from the field centers.) 


The last column of Table 2 contains additional notes where (ID) denotes possible 
identification with an object from another survey (see Table 3 and below); (S) suspect 
source which lies on an image processing artifact; (D) double source, either a blend of two 
comparable sources or two resolved peaks with a common envelope; (E) extended source, 
either an elongated source, or an object with an asymmetric or complex shape. In addition 
(D?) and (E?) indicate possible doubles and extended sources whose features appear at less 
than 8x o. Also included is the angular separation (arcminutes) between the peaks of double 
sources, the angular size of extended sources (measured from the peak to the most distant 
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feature, minus a slight resolution correction), and the position angle (degrees) measured from 
the brightest peak. We identified 124 double sources (45 of which have have both peaks listed 
individually in the catalog) and 174 extended objects with a range of sizes from 0.1’ — 2.2’, 
with 70% of the doubles and 60% of the extended sources being larger than 0.5’. 


Contour images of 54 source ‘s with interesting structure are included in Figure A4 and 
are indicated in the notes of Tavie 2. Positional corrections have not been applied to the 
maps and source positions are best taken from the catalog (Table 2). 


A4 Comparisons with Other Catalogs 


The VLA-NEP source catalog was compared with four previous radio catalogs which 
cover all or part of the VLA-NEP catalog area. These include the 1.4 GHz Green Bank 
catalog (White & Becker 1992), the 5 GHz Green Bank catalog (Becker et al. 1991) and the 
1.5 and 5 GHz VLA and Ryle Telescope observations of the 38 GHz 8C survey (Lacy et al. 
1992). A search radius of 160” and 50” was used for the 1.4 GHz and 5 GHz Greenbank 
surveys, respectively, equal to the uncertainty in those catalogs, and 30” was used for the 
8C catalog, which contains primarily lobe dominated objects. In addition we have made 
comparisons with the 2.7 GHz Effelsberg survey (Loiseau et al. 1988) using a variable search 
radius. For objects found by Loiseau et al. to be “pointlike” we used their stated error radius 
plus and additional 10”; for “slightly extended” sources the error radius plus an additional 
65", for objects with uncertain parameters we used 120” and for “extended” sources a search 
radius equal to the maximum source size Omax (arcseconds) was used, where 


1/2 


Omax = 10"(€ + 1) + 0.5 [(0, + A61)(82 + 402) — (4.35")?| (18) 


with ¢ the error code listed in Loiseau et al. (1992), 0, and 02 the fitted major and minor 
axes of elliptical gaussian fit to the sources, and A@, and AQ; their associated errors (all in 
arcseconds). 


VLA-NEP counterparts were found for all of the 1.4 GHz Green Bank sources and 
approximately 90% of the sources from the other three surveys. The comparison with the 
5 GHz Green Bank and 8C catalogs found relatively few with more than one VLA-NEP 
counterpart within the search radius (nine and three objects, respectively). However, 36% 
of the 2.7 GHz Effelsberg and 78% of the 1.4 GHz Greenbank sources have multiple VLA- 
NEP counterparts. This is partially due to the large error radius adopted for the latter two 
surveys, although many of these surveys have located sources which we have found to be 
doubles Jarge enough for the two lobes to be listed individually. Where the identification of 
a VLA-NEP source is not likely to be subject to confusion (> 80% of the flux arising from a 
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Figure A4 — Contour plots of 48 selected extended and double sources. The corrected peak 
flux density of the brightest peak on the map is given, as well as the lowest contour level as a 
percentage of the peak. Subsequent contour levels increase by factors of 2, with the highest 
contour level at 99%. One negative contour at the lowest level is also plotted. The maps are 
in order of increasing right ascension and all plots are at the same angular scale. The axes 
are not corrected for positional offsets. 
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Figure A4 — Continued. 
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Figure A4 — Continued. Contour plots of 6 sources with interesting structures using the 
same procedure as before except that the images are on a variety of angular scales. 


165 


single VLA-NEP object), we have determined spectral indices between the VLA-NEP flux 
densities and the 2.7 GHz and 5 GHz catalogs. Because of uncertanties in all the catalogs 
these identifications should cautiously and are only meant to serve as a guide. 


In addition to radio surveys, the VLA-NEP catalog was compared to catalogs of 
extragalactic objects listed in the NASA Extragalactic Database (NED). Using a search 
radius of 1’ for this correlation 42 galaxies (most brighter than 15 magnitudes) and 27 JRAS 
sources were found with VLA-NEP counterparts. VLA-NEP counterparts were also found 
for the galaxy cluster Abell 2280, the planetary nebula NGC 6543, and SN 1989P, although 
the latter is most likely a chance coincidence from the large search radius adopted. 


The results of these tentative source identifications are listed in Table 3. For each 
source listed in Table 2 with an ‘ID’ note the columns show: (1) the VLA-NEP source 
name (J2000); (2) the VLA-NEP total flux density; (3) the name of any 5 GHz counterpart 
from Becker et al. (1991) along with the cataloged flux density (mJy); (4) the name of any 
2.7 GHz Effelsberg counterpart from Loiseau et al. (1992) along with the cataloged flux 
density (mJy); (5) the name of any 1.4 GHz counterpart from White & Becker (1992) with 
the cataloged flux density (mJy); (6) the name of any 8C counterpart from Lacy et al. (1992); 
(7) the name of any NED counterpart that is not from the four previous radio surveys with a 
visual magnitude, if available from NED; (8) the spectral index between 1.5 GHz and 5 GHz; 
and (9) the spectral index between 1.5 GHz and 2.7 GHz; (10) a reference for sources with 
information from NED. 


A5 Source Counts 


A5.1 VLA-NEP Survey Results 


The VLA-NEP catalog gives an accurate determination of the sky density of 1.5 GHz 
sources in the flux density range 1 — 150 mJy, following the procedure described by Mitchell 
(1983) and Mitchell & Condon (1985). To minimize systematic errors due to uncertainties 
in the flux density correction parameters (see §A2), we have restricted the source counts | 
to objects within 15’ of the field centers. The two fields directly north of each mispointed 
field (F047 and F067) have also been excluded. We have included the region around the 
calibrator (F032), although the calibrator source itself (1749+701, VLA-NEP 1748.5+7005) 
has been excluded. Double sources with multiple catalog entries have been treated as single 
objects, with the peak being that of the brightest component. With these restrictions, we 
counted 1622 objects over a total sky area of 21.7 square degrees. Figure A5 shows the curve 
of 2( Spek), the integrated solid angle over which a source can be detected, vs. the peak flux 
density Speak in mJy. 
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Figure A5 — The effective sky area over which a source with a corrected peak flux density 
Sp ought to have been detected on our survey maps, as a function of Sp. 


The sky density of sources within a selected interval of total flux density is given by the 


number Nw of sources per square degree in that interval summed with weight [Q(Speax)]7', 
where Nw and its associated error ow are 


N 
D2 [2(Speat)]7* (19) 


Nw = 
oN 
(ow) = 2,12 Spear)] (20) 


The source counts are corrected to compensate for sources with total fluxes sufficient 
to have been included, but peak fluxes densities too low to have been detected and counted. 
To do this we have followed the procedure outlined by Mitchell & Condon (1985). The 
correction factor C; for sources in a given bin it is determined by considering a reference 
bin containing N brighter sources. The peak and total fluxes of the reference bin sources 
were scaled downwards to match the flux range within bin 1. The number Nniss of sources 
from the scaled reference bin that had peak fluxes below the 6 x o would not have been 
counted and were considered “missed”. If the reference bin contains all the sources within 
its flux range, then C; = N/(N — Nmiss). However, because the radio source population 
changes in the range of a few mJy, the reference bin chosen for bin i was the bin z + 1 with 
slightly higher fluxes (except for the lowest bin where the next two highest bins served as 
the reference). Because the reference bin itself may suffer from missed sources the correction 
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factor C; = NC,41/(N — Nmiss). Estimates of C were found for all bins with total fluxes 
below 10 mJy. The error associated with the correction o¢ = C//N, although where C = 1 
the associated error is undefined and we have not included it in the determination of the 
source count error. 


It is conventional to plot the source density distribution normalized to that expected 
in a static Euclidean universe: No = (S,'* — Sg**)/1.5, where Sp and Sy are the lower 
and upper boundaries of a selected flux density bin. The corrected, normalized differential 
number count and its error is then 


ns? = C (*) 2 (21) 
Cc 2 N 2 1/2 
Ons?25 = I() ory + (52) o| . (22) 


Bins are selected to include all sources with total flux densities greater than 0.3 mJy with 
successive bin boundries increased by factors of \/2. The lowest two bins were combined to 
compensate for the low number of sources and uncertanties in the flux determination. Finally, 
in order to better facilitate comparison with other surveys, the fluxes were scaled from 1.49 
to 1.40 GHz by assuming a spectral index a = —0.75. Table Al presents the resulting 
normalized differential 1.40 GHz counts, with the columns showing: (1) — (2) the boundaries 
of the flux density bins (mJy); (3) NV, the actual number of VLA-NEP sources found within 
each bin in the 21.6 squares degrees; (4) Nw, the weighted number of sources and ony 
the associated uncertainty; (5) C, the resolution correction factor and o¢ the associated 
uncertainty; (6) nS?°, the source count after a Euclidian normalization and ojs2s, the 
associated uncertainty. 


A5.2 The Radio Source Population 


The source counts in Table Al are plotted in Figure A6 with the flux densities plotted 
at the expected average flux 


= 1/7 


where 7 is defined as n(S) = kS~7 and we take y = 1.7 (Condon & Mitchell 1982; Windhorst, 
van Heerde, & Katgert 1984). 


For comparison, we have also included on Figure A6 other source count data obtained 
with the VLA (Condon & Mitchell 1984; and references therein), as well as the best fit curve 
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TABLE Al 
1.4 GHz Source Counts 


Flux Density 10-4 Nw n§25 
(mJy) (sr~? Jy***) 
a (4) (6) 
0.30 <S< 8.55 + 1.81 3.20 + 0.69 
0.60 <S< 8.48 + 1.20 6.47+ 1.03 
0.85 <S< 5.93 + 0.62 6.41 + 0.78 
1.2 <S< : 6.27 + 0.49 10.2 1.1 
7? <S< 2. 5.04 + 0.36 13.5 1.3 
24 <S< : 3.72 + 0.27 ate 16.1 1.2 
34 <S< : 2.75 +0.21 Sia 20.3 1.6 
48 <S< z 1.83 + 0.16 44 22.5 2.1 
<S< : 1.36 + 0.14 ore 28.4 3.0 
<S< 0.977 + 0.122 oes 33.9 4.2 
<S< 0.745 + 0.107 steiig 45.0 6.4 
<S< 0.634 + 0.098 ee 60.8 9.4 
<S< 0.529 + 0.089 tae 87.8 
<S< 0.438 + 0.081 ee 118 
<S< 0.317 + 0.069 ae 145 
<S< 0.181 + 0.052 ae 143 
<S< 0.166 + 0.050 he 222 
<S< 0.151 + 0.048 hes 341 
<S< 0.106 + 0.040 ies 395 
<S< 0.076 + 0.034 stave 476 
<S< 0.030 + 0.021 ets 319 
<S< 0.045 + 0.026 sei 805 


HoH OH HE HE HOHE HE HE EOE OE OE OE 


log (nS**) [sr* Jy*"*] 


~5 0 5 1 15 2 25 3 
log (Stos) [mJy] 


Figure A6 — The normalized differential source counts from the VLA-NEP survey 
extrapolated to 1.4 GHz. The VLA-NEP data (dark circles) are shown along with 
representative values from other VLA surveys (crosses) taken from Condon & Mitchell (1984), 
and the best fit curve derived from Westerbork surveys (Katgert, Oort, & Windhorst 1988). 


derived from a variety of Westerbork surveys (Katgert, Oort, & Windhorst 1988). We find 
good agreement between our results and these other studies, except at the lowest bin, where 
our source detection and flux determination is weakest, and for the higher flux bins, where 
the number of VLA-NEP sources is small. Within the range of 3 - 30 mJy we find our 
source counts falling slightly below the Westerbork line, although in agreement with some of 
the other VLA surveys. We do not expect significant instrumental effects in the VLA-NEP 
survey in this flux range, and suspect this slight dip may be a property of the radio source 
population in the direction of the NEP. A linear fit of the logarithmic VLA-NEP normalized 
differential source counts between 1 - 150 mJy gives a slope 0.68 + 0.03, consistent with the 
value of 0.71 + 0.03 derived by Windhorst et al. (1985) between 5 - 100 mJy. 


Several detailed studies have been made of the 1.4 GHz source count distribution (e.g., 
Condon 1984; Windhorst et al. 1985; Katgert et al. 1988; Windhorst et al. 1984 Windhorst, 
Mathis, & Neuschaefer 1990; Rowan-Robinson et al. 1993). These studies have shown that 
there is a shift in the nature of the radio source population near 1 mJy, with the brighter 
population being dominated by radio galaxies and quasars, and the fainter population 
by starbust galaxies. Therefore, the VLA-NEP survey should have found virtually all of 
the brighter population about the NEP, although only a fraction of the fainter, starburst 
population. 
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